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We study the spectral tilt of primordial perturbations in Hofava-Lifshitz cosmology. The uniform 
approximation, which is a generalization of the familiar Wentzel-Kramers-Brillouin (WKB) method, 
is employed to compute the spectral index both numerically and analytically in a closed form. We 
clarify how the spectral index depends on the inflation model and parameters in the modified 
dispersion relation. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

A power-counting renormalizable theory of gravity, 
proposed recently by Hofava [l|, has attracted much at- 
tention and generated a trend in quantum gravity. The 
essential aspect of the theory is broken Lorentz invari- 
ance in the ultraviolet (UV), where it exhibits a Lifshitz- 
like anisotropic scaling, t — > £ z t, x — > tx, with the dy- 
namical critical exponent 2 = 3. Modification to gen- 
eral relativity thus introduced in the high energy regime 
brings interesting consequences in cosmology [2], includ- 
ing the generation of chiral gravitational waves from in- 
flation 3], a bouncing scenario |4L a possible candidate 
for dark matter and others p. Since signals from 
the early universe could be observed for instance through 
the cosmic microwave background (CMB) temperature 
anisotropics and the stochastic gravitational wave back- 
ground, it is important to study the dynamics cosmo- 
logical perturbations in Hofava-Lifshitz (HL) cosmology 
(see, e.g., Ref. [7] and references therein). However, it 
should be noted that problems of HL gravity are still 
under debate [1,1,11. 

One of the basic ingredients in today's cosmology is the 
almost scale-invariant power spectrum of primordial per- 
turbations generated quantum mechanically in the early 
universe, and what is equally important is a small devi- 
ation from the exact scale- invariance. Inflation involves 
a natural mechanism to generate almost scale-invariant 
cosmological perturbations, and indeed this fact together 
with observations of the CMB and the large scale struc- 
ture (LSS) in the universe provides increasing support 
for inflation. The scale-invariance in this case essentially 
relies upon the almost constant Hubble rate, and a slight 
deviation from the exact de Sitter expansion is responsi- 
ble for the small tilt of the spectrum. 

Another interesting mechanism to generate scale- 
invariant cosmological perturbations was recently 
pointed out by Mukohyama |10| in connection with 
HL cosmology. (See also Refs. [ll OS 0, 03 03 El 
for cosmological perturbations in HL gravity.) The 
anisotropic scaling in the UV allows us to include higher 
spatial derivative terms in a scalar field Lagrangian, 
modifying the dispersion relation into the form uJ 2 cx fc 6 



in the early universe, which results in scale-invariant 
scalar field perturbations. In contrast to the standard 
mechanism, the spectrum in Mukohyama's scenario 
does not depend on the Hubble rate explicitly but 
is determined only through some model parameter. 
Therefore, the new mechanism liberates inflation models 
from those with a very flat potential. In Mukohyama's 
scenario, the dispersion relation which is not exactly of 
the form u> 2 cx k e will break the scale-invariance. 

In this paper we compute a slightly tilted spectrum of 
quantum fluctuations generated in HL cosmology, eval- 
uating accurately the s pec tral index by means of the 
uniform approximation jl7l . EEL El- The uniform ap- 
proximation is a generalized Wentzel-Kramers-Brillouin 
(WKB) method, which is mathematically controlled and 
systematically extendable. To the best of the authors' 
knowledge, the uniform approximation has not been used 
so far for models with modified dispersion relations. We 
consider a power-law/slow-roll inflationary background 
in HL cosmology, which provides a transparent example 
that allows for analytical and numerical calculations. Al- 
though our primary motivation comes from HL cosmol- 
ogy, we believe that the approach presented in this paper 
has a wide range of applications in different scenarios of 
cosmology. 



II. UNIFORM APPROXIMATION 

We consider a quantum field whose Fourier coefficients 
obey the equation of motion, 
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where the prime denotes the differentiation with respect 
to conformal time r], a(rf) is the scale factor, and (77) 
is the effective frequency 0, El- The field tfk will be 
a gravitational wave (multiplied by the scale factor) d 
and a scalar field perturbation (multiplied by a scale fac- 
tor) (Toj . Following the prescription of the uniform ap- 
proximation [I7I EH, El ; we rewrite the above equation 
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where we defined 



g(v) = -k 2 M + 
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and q(j]) = — 1/ 4r/ 2 , which guarantees the 
convergence of the approximation. This 
method provides a single approximate solution 
<Pk(v) = [ff(»?)/C(»?)]- 1/4 [a*Ai(C)+/3 fc Bi(C)], with 

2/3 



±(3/2)/; J±^n)dm 



where 



for 



C(r?) := ± 

rj > fj, — for rj < fj, respectively, and fj, which 
corresponds to the effective horizon-crossing time, 
is defined by g(fj) = 0, and Ai(z) and Bi(z) are 
the Airy functions. Choosing a k = \f^T^ arLC1 
Pk = — ia/7t/2, the mode function ifk is made 
to satisfy the Wronskian normalization condition 
ifk^fk' ~ VfcVfc = * an d has the asymptotic early-time 
behavior (p k (r]) — [2fc e ff(r?)]" 1/2 x exp[ 
for rj <C fj. 

Using the asymptotic formula of the Airy func- 
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tion Ai(z) 

(y / 7r) _1 z _1/4 e 22: "' /a for |z| — > oo with — 7r < argz < 7T, 
we find the asymptotic formula of the mode function after 
the horizon crossing: 
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Thus, we finally obtain the power spectrum P^ik^rj) = 

\ l Pk(v)/ a ( r ])\ 2 j an( i the spectral index of the quantum 
fluctuations can be estimated from dhiP^k, rj)/d\iik, 
which leads to 



lim 
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y/ g{m)dm- (5) 
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This provides a closed formula for the spectral index. 
Note, however, that this formula cannot be used for the 
cases in which the fluctuations reenter the horizon and 
the horizon crossing takes place again. 



III. POWER-LAW INFLATION 

For simplicity, let us consider power-law inflation [2d| , 
for which a(t) oc t 1+n in terms of proper time t and 
a(rj) = l/(— H*ri) 1+1 / n in terms of conformal time 
rj (—00 < rj < 0), with n > 0. We first study a sim- 
ple case, 



NPa? 



(6) 



In this case, it is easy to find the supcrhorizon behavior 
in the uniform approximation: 



i(r))Ae- v (3 + 2/n) l/ - 1 / 2 fc-( 1 +^ , (7) 



a := m 



I/-1/2 



1 + 3n/2 
/i + n(n + 1) ' 




FIG. 1: Growing mode solution to Eq. Q, obtained by em- 
ploying the uniform approximation, is compared to the exact 
result. The parameters are given by /i = 2, n — 3, H, — 1, 
and k 3 /M 2 = 10 -2 . The horizon-crossing time is 77 = —4.14. 



where A defines the typical amplitude of the fluctuation. 
In the de Sitter limit n — > 00 with the usual dispersion 
relation fj, = 0, we have .4 = For fi — 2 we find 

^4 = M, and thus the result of [10( is reproduced. For- 
tunately, for the dispersion relation ([6]), Eq. ([1]) can be 
solved exactly, and so we can compare the exact result 
with the approximate one. The appropriately normalized 
solution to Eq. ([TJ) is given by 
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(- ?7 ) 1 / 2 ij(l)(y( ?7 )), ( 8 ) 



where Hi, is the Hankel function of the first kind and 
y(yi) : = J v k e s(r)i)dr)i. This then implies that \tp k \ = 
C|<^oxact| at r] w 0, where C := V2^e- l/ i/ ,y - 1 / 2 /r(i/). In 
the de Sitter limit, n — > 00, for fi — we find C = 
3V2e" 3 / 2 ~ 0.95, while for fi = 2 we have C = ^ 
0.86 irrespective of n. It is important to note that C 
gives rise to a fc-independent correction to the amplitude, 
and therefore this does not affect the computation of the 
spectral index. The exact result and the approximate 
one are compared in Fig. [1] in the case of /i = 2. 

Now we are in position to study the spectral index for 
the above model. Eq. ([5]) or ([7]) gives 
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fj, + n(fi + 1) ' 



(9) 



which leads to the scale-invariant spectrum n ini = —3 
when /i = 2, as expected 10]. The result ((9]) is consistent 
with what was obtained in a different way by matching 
two asymptotic solutions at horizon crossing [l5j]. Since 
the slow-roll parameter e (:= —H'/aH 2 ) is explicitly 
given by e = 1/(1 + n) for power-law inflation, Eq. ([9|) 
can be written as n m i + 3 ~ — e(2 — /f)/(l + /i) to leading 
order in the slow-roll parameter. One can confirm that 
this reproduces the standard result for /1 = 0. Note also 
that n ini approaches —3 in the limit n — > 00 (e — > 0) irre- 
spective of the value /i. This means that in order to break 
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the scale invariance of the primordial spectrum we need 
the background universe which is not exactly de Sitter 
and /X 7^ 2. 

Next, let us study the case of HL cosmology, for which 
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In the case of gravitational waves, ot\ is non-zero Q, while 
a i = for a scalar field [ljj]- For this dispersion relation 
Eq. reduces to 



1 + n 



_ z -(2+n)/(l+n)^2 f ^y dz 
J (n- 1 + 3/2) 2 K -2 z -2n/(l+n) _ ' 



(11) 



where f(z) := 



Z 4 + Ct!Z 3 



a 2 z 2 + 1, 



(M/ J ff*)(fc/M) 1 /( 1+ ") and we defined a new variable z by 
-H*f] = (Mz/k) n /( 1+n \ The prime in Eq. JTTJ) stands 
for the differentiation with respect to z. The upper limit 
of integration z satisfies 



(1/n + 3/2) 2 K~ 2 z- 2n/{1+n) - f(z) = 0. 



(12) 



Figure [2] shows the contour of nj n i in the k and 1/n 
plane, obtained from a numerical integration of Eq. (jlip . 
The parameters are given by ai = and a 2 = 1 for the 
left panel, which corresponds to a scalar field perturba- 
tion of [10|, a>i = 2, a 2 = 1 and eci = —2, a 2 = E 
respectively, for the center and right panels. The latter 
two cases correspond to chiral gravitational waves ana- 
lyzed in 0. 

The result is understood in an analytic way as follows. 
One can recast the expression (fTTJj) in fc 2 ff oc fc 2+2 ^('J'' £ ) , 
where fi(r], k) — (l/2)dln//dlnz, and so (i is dependent 
on r) and k only through z. For the case k « 1, the 
solution to Eq. (|12[) is such that z ^> 1, and thus we 
havez- [K^^n- 1 + 3/2)]( 1+Tl )/( 2+3 ™) and n(k,f})~2- 
a l /{2z) + (a 2 - 2a 2 )/(2z 2 ). Then, Eq. © gives 
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2a 5 



2z 2 



1 



2 + 3n 



(13) 



This result applies to the case where the horizon crossing 
occurs in the UV regime, k 2 s ~ fc 6 . While, for the case 
k ^> 1, the solution to Eq. (fT5|) is given by z <C 1, which 
corresponds to the case where the horizon crossing takes 
place in the infrared (IR) regime, fc 2 ff ~ fc 2 . In this case, 
we may write z ~ [K _1 (n _1 + 3/2)]( 1+n )/" and fj,(k,ff) ~ 
a 2 z 2 + (3/2)aiz 3 . Then, Eq. © yields 



-3 1- 



a 2 z 2 + 3aiz 3 /2 



1 + 



The approximate formula reproduces the qualitative be- 
havior of the numerical result. The above results suggest 
that we must require small k (< 1) or large n (> 10), in 
order to have a spectrum that is almost scale-invariant 



but slightly deviated from the scale-invariant one. Note 
that 
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(1.5 x i057)V(i+«) Vl0 19 GeV 



M 



n/(l+n) 



10 15 GeV 



H* J VMpc" 1 



fc 



V(i+«) 



(15) 



The sign of ot\ signals different chiral modes of primor- 
dial gravitational waves [3j , and hence the spectral tilt of 
the gravitational waves is different for right-handed and 
left-handed modes. The above result shows explicitly the 
dependence of the spectral index on model parameters. 
Our result also gives a hint on the scenario of [l(|, in 
which quantum fluctuations of a scalar field can be the 
origin of curvature perturbations by promoting the for- 
mer to the latter, e.g., via the curvaton mechanism [21] . 



IV. SLOW-ROLL INFLATION 

Let us move on to more general slow-roll inflation mod- 
els. We continue to study Eq. ([1]), but now with a" /a 
replaced by a" /a -> C{i) 2 /r] 2 = 2a 2 H 2 (1 + E). Here, 
E is some general function whose explicit form depends 
on the problem under consideration: E = — e/2 for a ten- 
sor perturbations and a scalar field perturbation in the 
slow-roll regime; E ~ e — 3(5/2, where 5 (:= —(f>/H<f)) 
is another slow-roll parameter, for a scalar perturbation 
in conventional slow-roll inflation models, where <f> is the 
inflaton field and the dot denotes the differentiation with 
respect to the proper time. However, a further careful 
investigation is necessary for a scalar (curvature) pertur- 
bation in HL cosmology. For example, the fc 6 term in 
the dispersion relation disappears under the assumption 
of the detailed balance condition, according to the pre- 
scription of Ref. [13]. References @, [3 investigate a 
different case without the detailed balance condition but 
with the projectability condition. 

Assuming the dispersion relation ([5]), a similar calcu- 
lation for the right hand side of Eq. © yields n; n i + 3 ~ 
— (3e)/(l + n) + e/3 — 4E/3 to lowest order in the slow- 
roll parameters, where we assumed that E is of or- 
der of the slow-roll parameters and used the relation 
rj ~ l/[— Ha(l — e)]. In the case of n = 0, the well-known 
result is reproduced for the tensor and scalar perturba- 
tions. Similarly to the above, we can also derive formulas 
for the dispersion relation (I10p . Hereafter, we define z by 
the root of f(z) = (i? 2 /z 2 M 2 )(9/4 + IE- e/2), where H 
is the value of H at the horizon-crossing time. For the 
case M/H «: 1, we find z ~ (3F/2M) 1 / 3 > 1 and 

^ 1 + 3.-(l + ^ + -^)e + -e--£, (16) 

while for the case M/H > 1 we have z ~ (3H/2M) < 1 
and 

n- mi + 3 ~ -3 ( 1 - a 2 z 2 - ^z 3 ] e + ~e - -E. (17) 
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FIG. 2: (a: Left) Contours of niui on the 1/n and k plane. The parameters are ai = and Q2 = 1. The level of the contour 
is shown in the panel, (b: Center) Same as (a), but with ot\ = 2 and «2 = 1. (c: Right) Same as (a), but with ai = —2 and 
a% = 1. 



One can check that these results are consistent with those 
of the power-law inflationary background with large value 
of n + 1 (= 1/e), by taking the expression for the gravi- 
tational waves and the scalar field, i.e., E — — e/2. 

In Q, it was pointed out that the evolution of pri- 
mordial gravitational waves is different depending on the 
parity mode. From Eqs. (fl6|) and (fT7| . the difference of 
the spectral index is given by 



An ini = 



(18) 



»1) 



The difference might not be negligible for M/H -C 1, 
which corresponds to the case where the horizon crossing 
takes place in the UV regime (k 2 s ~ k 6 ). However, the 
difference is small for M/H ^> 1, which is the case where 
the horizon-crossing occurs in the IR regime [k 2 s ~ k 2 ). 
If we assume a scenario of Hofava-Lifshitz cosmology, 
one can predict the signal of gravitational waves, to be 
compared with observations (e.g., Ref. [Hj])- 

Now we consider curvature perturbations in the sce- 
nario of by setting ai = and E = —e/2. Provided 
that the spectral tilt of scalar field perturbations is di- 
rectly translated into that of curvature perturbations, we 
have 



27ea 2 H 2 , , - 

~M2 (M/H»l) 



(19) 



-2e + 



The CMB measurement by the WMAP satellite and LSS 
of galaxies provide a stringent constraint on the primor- 
dial spectral index. For example, the WMAP team re- 
ported n in i + 3 0.04 ± 0.013 based on the ACDM 



model [23j, combined with la supernovae observations. 
This gives a constraint on the inflation model and the 
parameters in the scalar field Lagrangian, assuming that 
the almost scale-invariant spectrum we observe is origi- 
nated from the Hofava-Lifshitz scalar field [lol ]. 



V. CONCLUSIONS 

In this paper, we have studied the spectral tilt of quan- 
tum fluctuations in HL cosmology, employing the uni- 
form approximation. In the power-law inflationary back- 
ground we have obtained the spectral index numerically 
and analytically. The deviation from the scale invariant 
spectrum is described by the parameters a% and 0L2 of 
the model, as well as the background evolution, i.e., n. 
The case of a general slow-roll inflationary background 
was also investigated. In order for Mukohyama's sce- 
nario 10] to be successful, the spectral index, given in 
terms of the model parameters by Eq. p^|) . is required 
to be 0.04. 
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